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CALCULATING MAX-EIGENVALUE AND MAX-EIGENVECTOR 
WITH JUMPS OF MATRICES 

ALI EBADIAN, SAEED HASHEMI SABABE, AND HOJR SHOKOH SALJOGHI 


Abstract. The eigenvector problem for an irreducible non negative matrix 
A = [aij] in the max algebras is the form A 0 x = Xx, where (A G x)i = 
max{aijXj), x = (xi, X2, ■ ■ ■, Xn)"^ and A refers to the maximum cycle geo¬ 
metric mean fJ,{A). In this paper, we exhibit a method to compute /r(A) and 
max-eigenvector by using mutations of matrices. Since the order of power 
method algorithm is 0(n®), the advantage of this paper present a faster pro¬ 
cedure. 


1. Introduction 

The max-algebra or (max, x) semi ring, is the set IR.+ , equipped with max as 
addition, and ordinary multiplication as its multiplication. It is traditional to use 
the notation © for max, and © for x . For vectors x = {xi),y = (yi) in R" and 
c € M+ the vectors x (By = (max{a:i, t/i}) and cx = (cxt) are defined element wise. 
The sum A © i? of two matrices is defined analogously. This structure satisfies all 
the semi ring axioms, i.e. © is associative, commutative, with zero element, © is 
associative, has a unit, distributes over ©, and zero is absorbing. This semi ring is 
commutative (a © 6 = 6 © a) , idempotent (a © a = a) , and non zero elements have 
an inverse for ©. The term diode is sometimes used for an idempotent semi ring. 
Using the new symbols © and © instead of the familiar max and x notation is the 
price to pay to easily handle all the familiar algebraic constructions. For instance, 
we will write, in the Max semi ring: 


■ 2 

0 ■ 

■ 10 ■ 


' 2© 10©0© 13 


■ 20 

II 

15 

13 


11© 10© 15© 13 


195 


(2) (3 + x)'^ = (3 + x)(3 + x) = 6 © 3x © = 6 © x^. 

Obviously with this multiplication for matrices, we can define max eigenvalues and 
eigenvectors of a matrix. If for some matrix A a vector x > 0 and a number y{A) 
exist such that 

(3) A©x = /i(A)x, (A © x)i = maxa^Xj, 

3 

Then x is called an max eigenvector of A, and /i(A) an max eigenvalue of A. For 
abbreviation we simply use eigenvalue and eigenvector instead of max eigenvalue 
and eigenvector. 

There are several applications of max algebras one can found in [niEiiniiii]. 
As it clearly showed in application, max algebras and especially eigenvalues and 
eigenvectors can describe some important properties. 

One of important application of max algebras is discrete events dynamical sys¬ 
tems (BEDS). In control engineering, a discrete event dynamic system (BEDS) is 
a discrete-state, event-driven system of which the state evolution depends entirely 
on the occurrence of asynchronous discrete events over time. Although similar 
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to continuous-variable dynamic systems (CVDS), BEDS consists solely of discrete 
state spaces and event-driven state transition mechanisms. We can correspond to 
the graph of every such systems, a matrix A which every entry Aij is the weight of 
arc (j,i) ID- 

A path of length I in a graph is a sequence of arcs 

( 4 ) 

A circuit of length I is a closed path of length I. A circuit which its first and last 
elements coincide is called a Cycle. 

Definition 1.1. A graph is called strongly connected if from any node i to any 
node j a path exists. If a graph is strongly connected, the corresponding matrix 
will be called irreducible. A non negative matrix A G M„ is said to be primitive if 
it is irreducible and has only one eigenvalue of maximum modules. 


Irreducible matrices have a special matrix form. Indeed, a reducible matrix 
A G Mn of size 1 is 0 and for n > 2, there is a perjump matrix P G Mn and an 
integer r with 1 < r < n — 1 such that 


( 5 ) 


P^AP 


B C 
0 D ’ 


Where B G Mr, D G M„_r, C G Mr^n-r and 0 G Mn-r,r is a zero matrix. A matrix 
A is irreducible if it is not reducible. 


Definition 1.2. The weight of a path —>■ I 2 —>■... I; —>■ ii+i is the sum of the 
weights of the individual arcs. The average weight of a path is its weight divided 

by the number of arcs: -|- H--|- — 1). The circuit mean is 

the average weight of a circuit. Any circuit of maximum average weight is called a 
critical circuit. Similarly, we can define circuit geometrical mean by 


( 6 ) 



X Ai3,i2 X 




Definition 1.3. A matrix A G Mn is said to have property SC if for every pair 
of distinct integers p, q with 1 < p,q < n there is a sequence of distinct integers 
fci = p,k 2 ,k 3 ,..., km-i, km = q, 1 < m < n, such that all of the matrix entries 
akik 2 ,ak 2 k 3 , ■ ■ ■ are non zero. 

Several related subjects about property SC can be found in m- We just mention 
to following theorem without proof. 


Theorem 1.4. Let A G Mn corresponding to graph r(A). The following are equiv¬ 
alent: 

(a) A is irreducible; 

(b) r(A) is strongly connected; 

(c) A has property SC. 


Braker and Olsder present in [1] a power algorithm to find max eigenvalue. They 
also found some useful properties. Esner and Van Den Driessche try to modify 
Braker’s method. In section 2, we try to present a method, based on concept of 
jump of matrices which is faster and has smaller absolute value of error. 

In section 3, we present a new method, based on max algebra, corresponding an 
optimization problem to a SR matrix such that minimize following function: 


( 7 ) 


e(w) = max 


a-ik - Wi/wk I 


^ik 


We also, introduce jumps of a matrix and its applications in symmetric, reciprocal 
and trasitive matrices and in the third section, we calculate the error refered to 
jumps and study behavior of jumps of transitive matrices. 
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2. Main Theorems 

To evaluate the determinant of a matrix A € Mn by Laplace expansion, we use 
the following formula: 

n n 

(8) det A = det Aij = ^(—det Aij, 

j=i i=i 

Where Aij is a matrix which obtain by eliminating the i—th row and j—th column 
of matrix A. This process is continued to evaluate det Aij while the size of remain 
matrix becomes equal to 2. To evaluate the determinant of a matrix with size 2, 
we simply differ the multiplication of entries on main diagonal and other diagonal. 
In fact in this process, we multiply on entry of every row and column. Determinant 
of the matrix is a special kind of linear combination of these multiplications. 

We can have an other formula to exhibit determinant of matrix A G M„. 

n 

(9) det A = '^sgn 

(7 i—1 

where the sum runs over all n\ perjumps a of the n items {1,..., n} and the ’’sign” 
or ’’signum” of a perjump cr, sgn a, is +1 or —1, according to whether the minimum 
number of transpositions, or pair wise interchanges, necessary to achieve it starting 
from {1,2,..., n} is even or odd. Thus, each product aicr(i)R 2 cr( 2 ) ■ • ■ ancr(n) enters 
into the determinant with a + sign if the perjump a is even or a — sign if it is odd. 

Definition 2.1. For matrix A G M„, every product of n entries of different rows 
and columns is called a Jump. A Jump with no main diagonal entry as its multipliers 
is called a Principal Jump. Other wise, we called it a Subordinate Jump. 

In this view, determinant of a matrix is a special linear combination of its Jumps. 

Example 2.2. Consider the following 3x3 matrix 

12 3' 

4 5 6. 

7 8 9 

We have 

(11) det A = lx5x9 — Ix6x8 + 2x4x9 — 2x6x7 + 3x4x8 — 3x5x7. 

In this example, 3x4x8 and 2x6x7 are principal Jumps and others are subordinate 
Jumps. 

Theorem 2.3. Let A be anxn non negative matrix. The part of jump that hasn’t 
any entry of principal diagonal form a cycle. 

Proof. Let a is a Jumps of A. Then entries of a are ..., which 

diagonal entries of A are eliminated. By definition, a is a cycle. □ 

Theorem 2.4. Let A be anxn non negative matrix. A part of a Jumps with non 
zero production of elements and without any element of a principal Jump form a 
cycle. 

Proof. It is not hard to show, using lemma 2. □ 

Definition 2.5. Let A be a square matrix. The set of all Jumps of A is denoted 
by M (A). Let a be a Jump of A. The product of non zero elements of a is denoted 
by p{a) and the number of non zero elements of a is denoted by Sa. 

Theorem 2.6. Let A be a non negative square matrix and a G M(A). Then 

(12) m(^) = max{ ^%/p{a) , a G M(A)}. 
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Proof. By using of two primer theorem and knowing that a circuit can appear on a 
vertex therefore element of principal diagonal can yield a circuit then whit attention 
to define of fi{A) proof is evident. □ 


Example 2.7. The eigenvalue problem for an irreducible non negative matrix 
A = [ttij] in the max algebras is 

(13) AiS)x = Xx (A'S) X)i = maxaijXj, 

3 

Where A denotes the maximum cycle mean fJ.{A). But by properties of max-algebras 


(14) 


J Uij 'ifihj) *5 in jump fJ,{A) 
\ 0 other wise 


And A*x = fj,{A)x we can simply obtain eigenvalues of A. 
Example 2.8. Consider matrix A as follows: 

0 8 1 ' 

3 0 2. 

4 1 1 



Then 


(16) det{A) = 0- 0- l- 0- 2- l- 8- 3- l-h8-2-4-hl-3-l-l-0-4, 

(17) p.{A) = max{ v'^, ^8 • 2 • 4, ^1 • 3 • 1} = 4. 

Where fJ,{A) is corresponding to subordinate jump 012 , 021 , 033 - Thus by equation 
A*x = fi{A)x we have 

■ 0 8 0 1 

(18) 0 0 2 [ xi X2 X 3 ] = [ xi X2 x^ Y' A. 

4 0 0 J 

So 

(19) [ Xi X2 X3 ] = [ 2 1 2 ] . 

Theorem 2.9. Let A be a n x n non negative matrix. Then 

(20) n an > max{p(a) , a G M{A)} /x(A) = max{aii}. 

Proof. When Iloii is greater than all p{a) ’s, then there one of the main diagonal 
entry has the greatest module. Then the maximum of geometric mean holds on the 
main diagonal. □ 

Theorem 2.10. Let A be a non negative sguare matrix. Then 
Va G M{A),p{a) < 1 4=^ h{A) < 1. 

Proof. It is not hard to see, since the maximum of geometric mean is equal to 
Sa — th root of some jump a. □ 

Theorem 2.11. Let A be a non negative sguare matrix. If maximum geometric 
mean fJ,{A) holds on a jump with only elements of principal diagonal, then max- 
eigenvector corresponding with fJ,{A) is Ci, where i is the index of fj,{A) = an, else 

X Ci. 
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Proof. By the condition of theorem, all entries of A* is equal to zero unless the 
entry corresponding to n{A) = an. Then by A*x = fi{A)x, Ci is the eigenvector 
corresponding to fJ.{A). If the maximum geometric mean fJ.{A) holds on any other 
jump, there is a non zero entry in every row of A* since there is an entry of every 
row and column of A in this jump. So A*x = fj,{A)x concludes fJ,{A) >0. □ 


The following theorem shows the relation between critical matrix and jumps. 


Theorem 2.12. Let A be a non negative matrix. The critical matrix of A is 
corresponding with a non zero jump which holds on its Jump cyclic geometric mean 
p.{A), therefore 


( 21 ) 


J Oij if (,Lj) exist in Jump of p,{A) 

\ 0 else 


Proof. It is not hard to see by [T^] and definition of critical matrix. 


□ 


3. Application of jumps in Producing Transitive Matrices 


An (entrywise) positive n x n matrix A = (a^) is called a symmetrically re¬ 
ciprocal matrix (SR- matrix) if aijOji = I for all i,j=I,...,n, thus an = 1. The 
firest time were introduced by Saaty, used in the analytic hierarchy process(AHP 
method) for multicriteria decision making. And An SR- matrix B = (6^) is 
called transitive if there is a positive n-vector W = {wi,W 2 ,... ,Wn) such that 

'Wi 

bij = — for I, j = I,..., n.thus sometimes it is required to deduce positive weights 

Wj 

wi,W 2 ,dots, Wn attached to the alternatives Ai,..., An respectively, from the SR- 
matix A. Approximating an SR-matrix A by a transitive matrix B is an importance 
step in the analytic hierarchy process(AHP) for decision making by attaching a 
ranking to the SR- matrix. In this way the alternatives can be ranked. For the 

Wi 

ideal case, an = — ior i, j = 1,... ,n. it should be noted a transitive matrix is SR- 
matrix of rank one. However, in a realistic case aij is only approximately given by 

'IVi 

—. one of the most important subjects in section is , costructing a weight vector 

Wj 

W = {wi,W 2 , ■. •, Wn) that there are several suggstions in literature for construct¬ 
ing a weight vector . Some of these methods is given below. Saaty proposes taking 
the perron vector of A . the vector W is chosen in such a way that the matrix with 
entries — has minimal distance from A in the Euclidean matix norm, i.e., that 


W'/' 

Wj 

is minimal, in the for costructing a weight vector W = {wi,W 2 , • ■ •, Wn) using that of 
max-eginvector of A that solves a useful optimization problem, namely minimizing 
the relative error useful optimization problem, namely minimizing relative error 


e{W) 


max 


aik - {wj/wj) 


^ik 


{*) 


but we in this peaper given new method for costructing a weight vector W = 
{wi,W 2 ,... ,Wn) with using of jumps of matrix . This method has the following 
advantages compared to other methods that are Short out the answer and Makes 
much more efficient error (*). the firest we mention several theorem about relation 
between jumps and SR-marix and transitive matrix. 


Theorem 3.1. A matrix A S M„(n > 3) is symmetrically reciprocal if and only if 
the production of all its subordinate jumps eguals 1. 
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Proof. Let A he a symmetrically reciprocal matrix. By definition, a subordinate 
jump contain some elements like aij,aji and an and by the properties of SR matrices 
on such a subordinate jump, = 1. Other side of the theorem is obvious. □ 

Theorem 3.2. A symmetrically reciprocal matrix A is transitive if and only if the 
production of all its jumps equals 1. 


Proof. Let Al be a transitive matrix, then there exist a wight vector w = {wi ,W 2 , ■. ■, Wn 
such that Oij = Wi/wj. Form of principal jumps is a^, Ojk, aki and form of subor¬ 
dinated jump is Oij, aji,aii. So we have 


( 22 ) 


Wi Wj Wk , 
- X X - = 1 


and 


Wi 


Wk 


Wi 


Other side of the theorem is obvious. 


w^ Wj Wi 

- X — X - = 1. 

Wj Wi Wi 


□ 


Lemma 3.3. For three integers a, &, c, following statements are equivalent: 

(a) < a6 < 1 -f c, 

(b) |a — 5| < ca and \^ — b\ < ^. 

Proof. See [12]. □ 

Theorem 3.4. Let A G Mn(JS.+ ) and we have its determinant by Laplace expansion. 
If Y\aij he the production of elements of each jump, w = {wi,W 2 , ■ • ■ ,Wn) be the 
weight vector and B = \bij] is a matrix such that bij = then following statements 
are equivalent: 

(a) T+c — ^ijbij ^ 1 c, 

(b) \aij — bij\ < coij and \ — — bij \ < —. 

In special case, if the maximum production is of a principal jump, then 

(23) c=yfp“-l. 

Where H ay w the maximum product of jumps and k is the number of elements of 
it. 

Proof. Let A be a SR matrix. By lemma [231 set a = aij and b = ^. So (a) and 
(b) are equivalent. 

But by © we have 

/r)^\ / \ I ^ik Wi/wk I yj „ 

(24) e('u;) = max - V w > 0, 

Wik 

And for every row of A and B we have 

'UUi 

(25) 1 + e(w) = uiaxtaik — )• 

Wk 

So 

(26) (1 -I- e{w))^ = (J^ max(aife —) = max(JJ Oij —). 

Wk Wk 

But max — = 1. Therefore 

Wk 

(27) (1 -I- e(w))^ = max (1 -I- e(w)) = ^ e{w) = — 1. 

□ 


Let A be a positive SR matrix and r > 0. We dehne = [alj] as follow: 

{ ai^j^T (Aoi) w the in the jump which ix{A) holds 

ajj^i^fr (ji;*i) w the in the jump which l//i(A) holds 

Oij other wise 
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Lemma 3.5. Let A be the matrix defined in [55] with n > 2. There exist positive 
numbers ti and T 2 such that /i(r) is strictly decreasing on (0,ti), is constant on 
iTi,T 2 ) and is strictly increasing on (t 2 , +oo). Moreover for t € (0,ti) U (r 2 , +c»), 
eigenvector is unique. 

Proof. Since /x(r) is constructed by elements of a principal jump and the product 
of principal jumps increases on interval (—oo,ti), so fJ-ir) decreases on this interval 
and increases on (t 2 ,oo). □ 

Lemma 3.6. By lemma 1531 if we define p,o = min{/i(r),r > 0} then following 
statements hold: 

(a) Ho = 1,T = 

(b) Ti < llaiHT2 > llaij, 

(c) r = 0, 

(d) r = To- 

Proof. (a) Since A is a SR matrix and by changing the elements of jump, the 
product of its element changes, but this minimum equals 1 if and only if A 
is transitive. So 

(29) /r(T) = = 1 ^ /r(r) = /jq = 1- 

(b) Let (io)Jo) be a coordinate which is an element of /i(T) = Vn Oij. 
Since by definition, changing the elements Oij to ^ and changing the coor¬ 
dinates to {jfi) cause change of elements to a^r, while Oij are in principal 
jump, so the maximum change holds when it occures on a principal jump. 

(c) It is clear by definitions. 

(d) If T = To Then A is a transitive matrix with = I. So rank [A) = 1. 

□ 

Theorem 3.7. Let ^(t) be the eigenvector defined on \2Hh nd element that changes 
aij, 0 < To < Ti and x,y be two max eigenvector of A{t) such that 

(30) A{to) ®x = h{t)x A{ti) iS)y = Hix)y. 

Then 

(31) h{to) < h{xi) => — > — i>2 and h{to) > nip) => — < — if^2. 

X\ Xi X2 Xi 

Proof. Let A is a SR matrix and ai ,2 is changed to ai^ 2 t than since To < ti and 
by attention of propety of max algebra and by helping of matix A(to)* and A(ti)* 
that introduced in second section we have 

Toai^2a;2 = HoXi 
^ik^k — 


^ps^s — 


and 


Tiai^2X2 = Hixi 
^ik^k — 


^ps^s — Ml^p 
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So we have the first equations 

Xi 

Toai,2X2 — /^O^l Mo — EoCi 2- 

X2 


and 


nai^2X2 = MiXi 


Xi 

fj,i = roai,2 — 

X2 


now since 

/ ^ yi 

fj,o < Ml ^ Toai,2— < roai,2 — 

X2 2/2 

X\ X\ 

therefor — < — and other Equations we can easily obtain Other inequality . 

X2 X2 

The second result follows from (i) by interchanging the roles of A{tq) and 


□ 


Corollary 3.8. With condition of Theorem VS.H j//i(To) = /i(Ti) then 


(32) 


Xi ~ X^ ~ X2 


Proof. It is a direct conclusion of previews theromes. 


□ 
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